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Abstract 

We suggest to compute structure functions in the Hamiltonian formalism on a momentum 
lattice using a physically motivated regularisation that links the total parton number to the 
lattice size. We show for the (f>\ theory that our method allows to describe continuum physics. 
The critical line and the renormalised mass spectrum close to that critical line are computed 
and scaling behaviour is observed in good agreement with the semi-analytical results of Liischer 
and Weisz and with other lattice simulations. We also demonstrate that our method is able 
to reproduce the Q 2 behaviour of deep inelastic structure functions and the typical peak at 
xb = 0. 

1 Introduction 

Hadron structure is probed by deep inelastic scattering (DIS). Much effort has been devoted by the- 
orists to compute quark or gluon distribution functions and proton structure functions from QCD 
with non perturbative methods. E.g., Martinelli et al.(l) have computed the first two moments 
of the pion structure function. These calculations are notoriously difficult; for the present status 
of lattice calculations based on the Euclidean action formulation see Ref. Q . In order to compute 
structure functions it seems natural to use the Hamiltonian formulation that allows to compute 
wave functions. E.g., scattering wave functions for glueball-like states in compact QED2+1 have 
been computed in a Hamiltonian formulation on a momentum lattice j3|. For a review of Hamilto- 
nian lattice methods seeQ[||. The usefulness of a momentum lattice to compute physics close to 
a critical point has been demonstrated in Ref. || Q| for the </> 4 -model. Unfortunately, Hamiltonian 
methods are known to pose problems when calculating the huge number of field degrees of freedom 
due to particle number. Nobody has yet succeeded in seeing scaling laws indicating continuum 
physics in a QC£>-Hamiltonian formulation, to the best of the authors' knowledge. In this work 
we shall demonstrate that these difficulties can be overcome. We introduce a new type of regulari- 
sation which excludes all parton momenta that we deem to be not directly related to deep inelastic 



1 



structure measurements. It turns out that such a regularisation scheme provides both an ultravio- 
let as well as a particle cutoff. We have chosen to test our method on the popular scalar (j) A theory 
in four space-time dimensions because this model has been thoroughly investigated numerically 
and analytically. The principal results are: 

1. Using a Hamiltonian approach we observe continuum physics (scaling laws close to the critical 
line) for the 4>\ and we see 

2. Altarelli-Parisi-type behaviour of </> 3 distribution functions in two, three and four dimensions 
and a peak at xb — similar to the peak occurring in quark distribution functions. 

2 DIS and a physical choice of Fock space 

Physical observables of relativistic field theories depend continuously on the experimental reso- 
lution. Every experiment has a finite resolution and coarse-grains the observed system. This 
corresponds cum grano salts to a regularisation of the field theory which describes the observed 
system. The hadron structure functions that we attempt to calculate are a prominent example. 
The most important experiment in order to probe the structure of hadrons is deep inelastic scat- 
tering (DIS) . Its simplest form is inclusive scattering of an unpolarised lepton off a hadronic target 
which we want to build our method upon. For the basics of DIS see, e.g.Q. Notation: The hadron 
in its ground state with four momentum P interacts with the probing lepton by the exchange of 
a virtual photon (our neutrino) with space-like four-momentum q. In Feynman's parton model, 
the impulse approximation is invoked. It is assumed that the proton consists of constituents, the 
partons, that are weakly bound when compared to the resolution ability Q := ^—q^q^ of the 
probing photon. In this approximation, the so-called Bj0rken scaling variable xb '■= 2 p ' q ^ can he- 
interpreted as the momentum fraction of the struck parton in the Breit frame (which is defined 
by the requirements that the photon energy q$ be zero and that the momentum q point in the 
opposite direction of the hadron momentum P). In this model, structure functions F(xb,Q) can 
be interpreted as some linear combinations of parton momentum distribution functions J{xb,Q)- 
For details see (8). 

We want to design a theory where just this picture emerges, describing only the coarse-grained 
particles that are probed in a DIS experiment. Consequently, we have to devise a regularisation 
that mimics the coarse-graining due to experiment as close as possible. In order to achieve this goal 
we have to decide on a conveniant Lorentz frame, because firstly the Hamiltonian is only defined 
with respect to such a frame and, secondly, a parton being a virtual particle is only defined with 
respect to the specific relativistic frame where experiment performs its coarse-graining of space. 
Seen from a different Lorentz frame, the parton would emerge as a jet of many partons because 
the boost generators of the Poincare group contain vertices that change the parton number! Our 
candidate is the Breit frame, of course. The struck parton in the Breit frame, however, is collinear 
to the hadron momentum. Thus, we have to allow for non-collinear partons in more than one space 
dimension. Therefore we use a 'generalised Breit frame' where q — but q does not have to be 
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collinear to P. Please note that only in such a 'generalised Breit frame' is Q = \q\ and only there 
can the wave length A = of the virtual photon be interpreted as the resolution ability of the 
experiment! Please note also that the momentum P strongly restricts the domain of all possible 
photon momenta q in the generalised Breit frame. We require therefore that the 'observable ' Fock 
space consist only of partons whose momenta can (but need not) reverse momentum — p = p + q 
if struck by a virtual photon in the generalised Breit frame. Parton momenta that can not be 
related by p — —q/2 to some virtual photon with momentum q are skipped. Then all 'observable' 
momenta of the ground state lie inside a sphere 

(p-P/2f<\P/2\ 2 (1) 

~* ~* O 2 O 2 11 

around P/2 with diameter \P\ because xb — jp^ — ~ p~q ^ es between and 1. 
If we go onto a momentum space lattice with momenta k := nAk where ft is an integer vector and 
Afc is the momentum lattice spacing, then the constraint (Q) on the parton momenta provides an 
ultraviolet as well as a particle cutoff. This holds because only a finite number of partons with 
positive longitudinal momentum p e := p ■ e can share the positive total longitudinal momentum 
P ■ e, where e := P/\P\ denotes the direction of the total hadron. If one performs the continuum 
limit, the system's velocity v = P/E grows and one ends up with the velocity of light. For a 
parton with momentum fraction xb = const the continuum limit a — > also implies Q = ^2 — -> oo 
and hence is equivalent to the Bj0rkcn limit. By the way, the transverse momenta in this limit are 
suppressed by the increasingly dominant kinetic term, at least in asymptotically free and trivial 
theories. 

Since we want to calculate the critical line - which depends on the regularisation - we had to 
ensure that our regularisation([l|) be well comparable to to the lattice regularisation of ref.j9j. 
To this aim we started on the usual lattice with the position space lattice spacing a — N ^ k as 
ultraviolet regulator, where 2N + 1 denotes the number of lattice points in each space direction, 
and put the total momentum P on one of the edges P — (1, 1, l)AkN of the momentum lattice. 
If we then apply our regularisation(|l|), every component of the parton momenta is positive, i.e the 
combination of our regularisation with the lattice cutoff restricts the parton momenta to the first 
octant of the momentum lattice which further simplifies the set-up of the fock space. The maximal 
number of 'partons' equals N in this description. We ought to add that for a small lattice such as 
2N +1 = 7 the additional lattice spacing cutoff does not exclude states that are allowed by our 
regularisation - so this seems to be o.k.. 

A further feature of the proposed ansatz is remarkable: the connection between boosts and renor- 
malisation group transformations. A boost that reduces the total momentum \P\ equally reduces 
the number of observable Fock states in our approach and hence constitutes a coarse-graining 
procedure on a momentum grid with fixed infrared regularisation- in contrast to that, a boost 
without our regularisation would increase the (unobservable) particle number. Since the hadron 
momentum \P\ = — is the cutoff in our approach, the scale 1/a is reduced and the bare parameters 
have to be re-adjusted. This means that such a boost acts in a way similar to a renormalisation 
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group transformation in the Kadanoff/ Wilson sense. Such a transformation is not invertible. In- 
tuitively this corresponds to the fact that decreasing the momentum of the probing lepton relative 
to the target hadron corresponds to a lesser structure resolution. If \P\ gets too small we have 
coarse-grained too much and the system breaks down, of course. In the momentum-sector with 
P = only one state, the vacuum, remains. Within our method we cannot describe a hadron 
in its rest frame, just as one cannot perform structure measurements without a particle beam 
moving relative to the target. This means that we cannot give a Fock space expansion of the 
physical vacuum either. This is no draw-back, it is a desirable feature. It makes simply no sense 
to calculate structure functions for the vacuum. We treat it as a background field with the same 
quantum numbers as the perturbative vacuum because only vacuum expectation values of fields or 
of the energy (Casimir effect) originate from the vacuum. The fact that the perturbative vacuum 
is annihilated by our Hamiltonian does not mean that it is the physical vacuum -the latter one 
is far from being trivial. It only means that per constructionem our theory is designed for deep 
inelastic scattering and breaks down if we set P equal zero. If we, however, boost the momentum 
grid also, the transformation is an invertible scale transformation making the physical quantities 
change according to the Callan-Symanzik equation. Three little remarks at the end of this section: 

1. Theoretically, our regularisation(|l|) admits partons with zero momentum. An infinite number 
of these can be present because they do not take away longitudinal momentum. It is these that are 
responsible for vacuum expectation values. Fortunately, for a system with spontaneous symmetry 
breaking as the <f> theory, we can treat them on the classical level since at least in the infinite 
volume limit, quantum mechanical fluctuations of the vacuum expectation values must vanish. To 
describe the axial anomaly of a non-abilian gauge theories, however, we shall have to allow for 
dynamic gauge fields with zero momentum in order to describe fluctuations of the Dirac sea. 

2. It is obvious that our approach sheds some light on the relation of the so-called front- and 
instant forms|nj in two dimensions. A discussion of this must be reserved for another paper. 

3. Our method would, of course, permit to calculate the hadronic tensor and thus the structure 
functions directly from the electromagnetic courrent operators. But the calculation of distribution 
functions is much easier. 

3 Application on the scalar (f)\ theory 

Now we apply our method to the (f>\ theory as a testing ground. We shall compare it to some 
existing results |||| in the lattice regularisation. The Hamiltonian of the 4 theory is 

H=\J d 3 x(l/2((f ) 2 + (^) 2 + mW) + f 4 ) (2) 

where 4>{x) is the scalar field, mo and go are the bare mass and coupling constant respectively. 
We set the scale Ak to one and insert the discretised and quantised solution 

<j>{x) = 1 =(a{k)e- lk ^ + a ] {k)e +ik ^) (3) 

S ^fj2uj{k) 
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of the kinetic field equations of motion back into the Hamiltonian, where a, are the usual 
annihilation and creation operators [a^(fc),a(Z)] = 8% f respectively and oj(k) — ko — \J m\ + k 2 
is the kinetic energy of a parton. This yields the discretised Hamiltonian that we diagonalised. 
We denote the n-th eigenvalue of the Hamiltonian Matrix as E n and the n-th physical mass as 
M n := \J E 2 — P 2 . We did not normal order - which would not have been allowed for an interacting 
Hamiltonian - but we subtracted the vacuum energy. The latter is equivalent to the requirement 
that the perturbative vacuum be annihilated by the Hamiltonian. Only if the vacuum-energy 
Eq is gauged to be zero can we use the mass-shell condition M n := *J E 2 — P 2 to calculate the 
physical mass spectrum from the energy spectrum. We do not have to take the field renormalisation 
constant into account, but one problem arises. The critical line lies entirely in the region where 
to 2 , is negative and we cannot build up the Fock-space in terms of partons with imaginary masses. 
Instead, we have chosen to build it up in terms of partons with mass mjf by splitting the TOq0 2 /2 
mass term into a kinetic m 2 K 4> 2 /2 and a dynamic m 2 nt (j) 2 /2 part where m 2 K = m 2 , — m 2 nt > 
is the kinetic energy of the partons and m 2 nt is a sufficiently large negative number that takes 
the role of a coupling constant. We have presented our results for the choice m 2 K = 0. We have 
found, however, that the critical behaviour was almost independent of a specific choice of in 
the N = 4-sector, if only m 2 K is reasonably small. 



4 Critical behaviour of the scalar 4 -theory 

We diagonalised the described Hamiltonian on a 7 3 and a 9 3 lattice. To do this, Fock spaces of 
only 6 resp. 21 states are needed(and even less if the symmetry is removed in the symmetric 
phase). In order to compare our results to the work of Liischer and Weiszj^] and to the work of 
D. Berube et al. |J we have to introduce the parameters A and k which are related to too and 
A by m\ = (1 — 2X)/k — 8 and go = 6^y. In Fig.l the first of the three data columns of the 
renormalised masses run in terms of K as given by Liischer and Weisz is shown and compared to 
our ground state masses M\. We have diagonalised the Hamiltonian for 50 x 30 different bare 
parameters mo, go and fitted the results to obtain the energy spectra in terms of A and n. It can 
be seen in Fig. 1 that the shape of the function is well described even on the 7 3 -lattice. For the 
9 3 -lattice we are able to describe most of the scaling region Mi < The ground state mass 
M\(k) inside the validity domain l/(2N) << aMi « 1 of our model is seen to obey the scaling 
law M\ ~ C , 4T 1 / 2 |Znr|~ 1 / 6 , where C4 is a constant^ and r := 1 — n/n cr it denotes the reduced 
temperature, meaning that we have correctly described the critical exponent v — 1/2 plus the 
logarithmic correction |?nr| -1 / 6 in our Hamiltonian approach. It can also be seen that the critical 
point is well reproduced. This is true also for curves of different A close to the Gaussian fixed 
point as it is shown in Table 1. where we listed the quotients ajy := K C ritl k Jhu f ne quantity 
denotes the point where the renormalised mass Mi becomes imaginary(It becomes real again in the 
asymmetric phase if we break the Zi symmetry by choosing a non zero field expectation value.). 
Table 1 shows that the quantity n^} t is an excellent estimator for the critical line. In Fig. 2 finally, 
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the mass spectrum related to the first excited state, M n /Mi, is depicted for fixed A in terms of k. 
It clearly demonstrates the fact that the partons of the <f>\ theory interact by repulsive forces. The 
first state above the vacuum is a one particle state with negative parity and the positive parity 
state above consists of two particles. 

5 Distribution functions 

Since the <j> 4 theory describes partons with repulsive two particle exchange forces and since it 
has a trivial, i.e. free continuum limit, its distribution functions are not comparable to the QCD 
structure functions. In the 4 theory, a renormalisation group transformation(RGT) decreases the 
lattice size as well as the marginal coupling constant. In an asymptotically free theory such as 
QCD, the coupling constant increases under such a transformation. This means a RGT enforces 
particle creation for QCD, whereas particle creation is weakened for the scalar theory. Since 
the renormalisation group in the Wilson/Kadanoff sense is used to make a complicated system 
numerically tractable by coarse-graining it into a theory with equivalent long distance behaviour, 
this means for the </> 4 -theory that the coarse-grained system has a strongly reduced particle content 
so that a particle cutoff may be applied after a few RGTs. This can not be expected to hold for QCD 
where a RGT actually increases the interaction. This effect slows down or may even invert the net 
particle number decimation of a RGT due to a reduction of momentum space points and this is the 
problem we designed our method for. The distribution function f(x B ,Q) of finding some parton 
with momentum fraction xb inside the hadron is determined by the parton momentum distribution 
function f(p, P) for finding a parton with momentum p inside the hadron with momentum P. 

Since Q is a dimensionful quantity, its scale is set by the lattice spacing Q ~ — and thus 

a(m ,go) 

depends on the bare parameters if one keeps the renormaliscd mass fixed. The continuum limit 
a — > corresponds to the the limit towards arbitrarily high resolution ability. If one keeps the 
rcnormalised mass and the renormaliscd coupling constant fixed, then Q is a function of the bare 
coupling constant go and vice versa - invertibility of Q(go) assumed. Hence f(xs,Q) is related 
by the function Q(go) to the distribution function f{xB,9o) that only depends on dimcnsionless 
parameters. Consequently, a calculation of the distribution function along a renormalisation group 
trajectory can be used to compute some peculiarities of the Q-dependence of the quark structure 
functions such as the increasing of the peak at ib =0 when the resolution Q is increased. In 
order to demonstrate that the suggested method is capable of describing them, we calculate the 
distribution functions of the (j) 3 theory whose forces are attractive one particle exchange forces. 
Although this theory is known to suffer from an unstable vacuum, it serves well for our purpose of 
demonstrating qualitatively that the singularity at xb =0 and its change in terms of the resolution 
Q of experimentally measured structure functions can be modelled in our approach (The unstable 
vacuum of the 3 theory does not allow us to calculate meaningful ground state masses that are 
needed to specify renormalisation group trajectories and consequently the exact relation between 
the resolution Q and the bare coupling constant g . But the distribution function J{xb, go) suffices 
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for our purposes). To increase the number of particles that can be modeled on our tiny workstation, 
we have chosen to present the results of the two dimensional 4> 3 theory. The results, however, are 
qualitatively alike in higher dimensions except for the smaller particle number. In Fig. 3 we depict 
the function f(xB,go) to illustrate that our method is - at least in principle - powerful enough to 
describe the known features of the proton structure function. This can by no means be taken for 
granted. Fock space methods that have to apply a low particle number cutoff in addition to the 
phase space cutoff are in principle unable to show this feature |jl(J. This is so because a system of 
n identical observable particles must have an expectation value of xb around < xb >= l/n for 
symmetry reasons. Furthermore, the behaviour of the structure functions in terms of the resolution 
Q is dominated not by the power law rest effects of the interaction between the partons but rather 
by the fact that a higher resolution power unveils more particles to the observer. With a low 
particle number cutoff, the splitting functions appearing in the Altarelli-Parisi equations would be 
meaningless, anyway. For the (j) 4 theory, however, an additional particle cutoff is perfectly feasible 
for describing the lowest-lying eigenstates of the Hamiltonian. 

6 Discussion 

We have shown that the suggested principle that only measurable quantities be described works 
well for the scalar theory. Even a 6 x 6 Hamiltonian matrix is able to correctly describe important 
features of the critical behaviour of the </) 4 -theory. We are, however, aware of the fact that there is 
a huge difference between the scalar </>|-theory which does not even possess a nontrivial continuum 
limit and nonabilian gauge theories which do possess a continuum limit with strongly bound 
particles and also an additional topological structure. But we have reasons for being optimistic. 
Our method has been designed to sooth one major nuisance of highly relativistic systems, the 
uncontrolled particle production. This feature has not even been necessary for the (p 4 theory where 
we could introduce a modest particle cutoff without mutilating the results. This is encouraging 
since we did not aim at the scalar theory in the first place. The aim is the calculation of structure 
functions, one of the outstanding problems of particle physics, directly from QCD where high 
particle numbers are inevitable. We hope that our results can be extended to gauge theories 
and stress that the suggested method has the power to significantly simplify the Kogut Susskind 
Hamiltonian because it also works on a position space lattice. This is important since there is no 
QCD Hamiltonian for momentum space available yet. The Kogut-Susskind Hamiltonian being a 
position space method is plagued by mirror fermions. Our approach, taking only parton momenta 
of the first octant in the brillouin zone, formally eliminates them from the set of 'observable' 
partons, which seems to be a remedy to this problem (We break the mirror-symmetry by the 
choice of a specific hadron momentum P different from 0). If this suffices to describe the axial 
anomaly remains to be seen. We have to be very cautious here: the Nielsen-Ninomiya no-go 
theorem [n2| has yet always resisted attempts to circumnavigate it. 
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Figure Caption 

Fig.l The ground state mass tor(A, k) in lattice units (a = 1) is depicted as a function of the 
lattice parameter k for the the semi-analytical data of ref. 0] (points) and for our data 
obtained on the 7 3 -lattice(dashed line) and on the 9 3 -lattice(solid line). The parameter A is 
taken to be A = 0.00345739. 

Fig. 2 The lowest lying mass spectrum on the 9 3 lattice is depicted in terms of k for A = 0.00345739. 
The ground state is set to one. 

Fig. 3 The distribution function f(xB,go) of (/>?> in terms of the momentum fraction Xb and the 
coupling constant g$. The mass mo is set to be Too = 3Afc. The maximum particle number 
N is N = 11. 
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8 Table Caption 



A 


.001 


.005 


.01 


.02 


.03 


.04 


.07 




0.125202 


0.125991 


0.126968 


0.128604 


0.130096 


0.133096 


0.133825 


a 3 


1.01327 


1.01194 


1.01192 


1.01718 


1.0227 


1.02664 


1.03075 




1.00949 


1.00518 


1.0043 


1.00262 


1.00079 


1.00068 


0.95589 



The semi-analytically determined critical points K cr i t in terms of the lattice parameter A and the 
quotients an ■= Kcrit / ' K cr2- The quantity K^t denotes the points we found where the renormalised 
mass M\ becomes imaginary. A tiny systematical error of about 1% can be expected to be present 
due to effects of our having interpolated the function Mi (mo, go) on a finite (but large) number of 
parameter points. 
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